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ABSTRACT 



The use of an electronic analog computer in solving the differential 
equation of motion of the free piston gas generator has been investigated. 
To accomplish this end, the free piston has been treated as a specialized 
mass-spring system, wherein the forces were derived from thermodynamic 
considerations rather then from elastic considerations. The first law of 
thermodynamics was solved on the computer three concurrent times in order 
to produce electronic analogs of an internal combustion engine, a recipro- 
cating compressor, and a gas spring. Coulomb and viscous friction were 
also simulated. The results of these four analogs were then combined and 
integrated as indicated by Newton's second law. The computer yielded a 
solution which represented the piston motion,, A total of 26 amplifiers, 
seven diode limiters, and six electronic multipliers were used to solve 
the problem,, 

The characteristics of an existing free piston machine were intro- 
duced into the computer as a final check of the analog. The accuracy of 
some of the results was less than desired due to errors in a two-quadrant 
electronic multiplication involved in the gas spring analog. The other 
three analogs performed very satisfactorily. For accurate and repro- 
ducible results, some other method, such as direct function generation 
or single quadrant multiplication, must be used in the simulation of 
the gas spring. 
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1. Introduction 



A great deal of effort has been put into the solution of the free 
piston engine problem. All of the methods currently in use are digital 
methods involving a point-by-point calculation of the variables asso- y 
dated with such a device 0 These methods are time-consuming and cumber- 
some due to both the trial and error element and the many geometric and 
operating parameters involved. In addition, the digital methods can be 
made trao table only by incorporating many numerical assumptions! some 
of which introduce considerable error. 

London^ ^ has developed a " Thermodynamic-dynamic analysis" whereby 
the calculated parameters of a "standard" machine may be scaled by 
certain affinity relationships to establish the design of larger or 
smaller machines. This method eases somewhat the numorioal burden of 
calculation, but the results proved to be quite sensitive to the assump- 
tions made; and in order to predict the characteristics of an existing 
machine, the assumptions were relaxed in four successive groups before 
the desired accuracy was obtained. This procedure still presented a 
formidable numerical task. 

In addition, a solution by digital methods becomes more complex in 
the investigation of part-load characteristics. 

The purpose of this investigation is to determine the possibility 
of solving the differential equation of motion of the free piston on 
the Boeing Electronic Analog Computer (BEAC), whioh is described in 
Appendix VI. If this oan be done, the oomputer will electronically 
simulate the free piston, presenting as a solution the piston motion, 
with most of the pertinent quantities also available. The advantages 
of this method include j 

1. The feedback principle, typical of analog computation. 
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eliminates trial and error. 

l 

2. All of the parameters may be taken into consideration, and 
the effect of each may be observed. 

3. The solution is immediate, thus permitting the study of many 
combinations of parameters in a short time. 

4* Fewer simplifying assumptions are required, and these may be 
relaxed to any degree compatible with the size of the computer* 

5* Part-load behavior can be studied as easily as full-load 

behavior* 

6* Time may be scaled to permit s low-motion studies* 
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2* General Background of Free Piston Systems* 

Figure 1 illustrates schematically two of the possible free piston 
configurations. Figure 1(a) illustrates the free piston air compressor, 
where the useful output is compressed air; and Fig, 1(b) illustrates the 
free piston gas generator, where the useful output is high pressure- 
temperature gas for use in a turbine* 

This investigation will encompass only the free piston gas generator, 
A slightly larger computer is required to simulate the free piston com- 
pressor* 

The free piston machine depicted in Fig* 1(b) consists of two 
opposed pistons, the inboard ends of which operate as a two stroke diesel 
engine* Another pant of each piston acts as a reciprocating air com- 
pressor* As the pistons move outward, the expansion work of the engine 
and compressor is stored as energy in the bounce chambers* This energy 
is then used on the inward stroke for the compression and delivery of 
the air in the compressor and the compression of the engine charge* 

The outstanding features of this free piston machine arei 

1* Simplicity of construction due to absence of crank 
mechanisms* 

2* Opposed piston configuration eliminates vibration* 

3* No piston side thrust arising from crank inertia forces* 

In the case of the gas generator, it would be well to review the 
features of a typical cycle* 

Referring to Fig* 1(b), which illustrates, say, the pistons just 
starting the instroke. Beginning at this point, the following successive 
processes may be visualized: 

1* High pressure in the bounce space forcing the pistons 
together; compressing and delivering the air in the compressor space, and 
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compressing the engine charge when the scavenge ports are oovered. 

2. Reversal of piston motion occurring during injection and 
burning of fuel, 

3, The expansion of the fluid in the engine cylinder and the 
intake stroke of the compressor; while the bounce cylinder is storing the 
expansion energy of the engine and compressor. 

4. Opening of the engine ports during the latter portion of 
the outward stroke, delivering working fluid to turbine, 

5, Pressure in bounce space causes reversal of piston motion, 
and cycle repeats. 

The preceding cycle is represented on Temperature -Entropy coordi- 
nates in Fig. 2, It is to be noted that the processes involve different 
masses, and can be represented only on a '’per pound'* basis. This diagram 
is similar to that in common use in free piston literature. It should be 
pointed out that the following two assumptions are made in producing this 
diagram; 

1. The scavenge air is furnished at oonstant pressure and 
temperature, related to the compressor d is charge state. 

2. Delivery to the turbine occurs at constant pressure and 
terape ra ture , 

These two assumptions are intimately connected with the mechanioal design 
of the engine ports and the physical size of the turbine receiver. It 
is beyond the scope of this investigation to oonsider such problems. The 
foregoing, therefore, constitute the first two assumptions in this 
investigation. It must be pointed out, however, that these assumptions 
must be relaxed if it is desired to investigate transient behavior, inas- 
much as these particular parameters are the primary transients. Before 
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Tempera tar* 




Process 1-2 - Compression in Compressor 0 
Process 2-3 - Pressure Drop into Engine 0 
Process 3-4 - Compression in Engine 0 
Process 4-5 - Burning of Fuelo 
Process 5-6 - Expansion in Engine 0 
Process 6-7 - Blowdown as Exhaust Ports 0pen o 
Process 3»7«=8 - Mixing and Scavenging, 

Process 8-9 - Expansion in Turbine 0 



Figure 2, Temperature - Entropy Representation of Fres Pistcn 
Generator- Turbine System, 
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proceeding further, the following are some additional assumptions* 
covering the device as a whole: 

1. That the compressor intake is at 70°F. and 14 # 7 psia. 

2. That the operating fluid obeys the ideal gas laws# 

3. A nominal pressure drop of five per cent across each valve 
and engine port. This is arbitrary, and may be refined to any desired 
value . 

4. That the geometry is symmetrical with respect to a plane 
passing through the center of the machine perpendicular to the hori- 
zontal axis. The end projection of this plane will be referred to as 
the "oenter line." 



t 



i 



♦Refer to Appendix I for all of the assumptions used. 
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3 0 Tha General analog of the Free Piston 



From a dynamic point of view, the piston in a free piston machine may 
be considered as a rather special sort of mass-spring arrangement, wherein 
the "springs" yield forces derived from thermodynamic considerations 
rather than from elastic considerations. 




Figure 3, Spring Analogy of the Free Piston, 

The system may be visualized as similar to that depicted in Fig, 3 above, 
where the spring forces are proportional to the conventional Pressure- 
Volume representation of the appropriate space. Not shown, but also 
present, is the force of static and viscous friction. Qualitative 
representations of these forces appear in Fig, 4. 

To illustrate the analog computer approach to this tyoe of problem, 
consider the simple mass-spring arrangement illustrated in Fig, 5. The 
differential equation of this motion may be immediately written down as 
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Volume or Displacement 

(c<>) Bounce Cylinder (d 0 ) Friction 

Figure 4o Typical Forces Acting on the Free Pis ton 0 




Figure 5. Simple Spring-Mass System# 
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771? + Kf + V f - O ) 



where m ia the mass 

v is the damping oonstant 
k is the spring oonstant 
i is the acceleration 
j is the velocity 
jP is the displacement 



By solving for the highest 
derivative. 




Recalling that the operational amplifier can add and/or integrate 
with respect to time, a computer soheraatic may be immediately devised to 
solve this problem. This is illustrated below in Fig. 6. 




Figure 6. Electronic Analog of the Simple Spring-Mass System. 



^Refer to Appendix II for analog computer table of symbols. 
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This simple problem points out the utility of analog feedback. That is 
to say, the immediate stabilization to the correct solution by using this 
solution to generate the derivatives, which are then integrated to form 
the solution again. 

Turning now to the free piston problem, and referring to Fig. 3; the 
equation for this motion may be formulated in accordance with Newton’s 
second law as 

MX = Pg(X,X) + F 0 (X,X) - FbU,X) " f F( X ), 3d 

whore M is the piston mass 

'xX ) is the force on the engine piston 

Fo(z^X ) is the force on the compressor piston 

F^^jX) is the force on the bounce piston 

F^(x) is the force due to friotion 

x is the piston displacement 

X is the piston velocity 

X is the piston acceleration 

In a manner similar to the simple mass-spring problem, a computing ar- 
rangement can be immediately proposed. Suoh an arrangement is outlined 
in Fig. 7. In this figure the blook representing the engine, compressor, 
bounce, and friction, are rather complex computing networks with many 
more inputs than shown. The overall principle, however, is basically 
the same as that outlined in Fig. 6. 

This fact has two very important consequences in that : 

1. A completely compatible set of parameters must exist in 
order for the computer to produce a sustained oscillation. If this con- 
dition is not met, a situation exists which is analagous to overdamping 
or negative damping of the simple system. 

2. Sinoe equation 3,1 also expresses the work-energy equality, 
a solution of this equation will automatically satisfy the condition 
that the engine work appears as the sum of compressor and friotion work 0 

The computer arrangement as shown in Fig, 7 will produoe displace- 
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Figure 7. Basic Computer Arrangement for the free Piston 
Problem,, 
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ment as a function of time. Also available will be the various foroes 



or pressures contributing to the motion. Related items such as air-fuel 
ratio, scavenge ratio, and power output will not be available direotly 
from the computer but oan be easily caloulated. 

The inputs to the system will be : 

1. Geometrio parameters, which will be fixed for a given 

maohine . 

2. Initial pressure conditions, which arei 

a. Intake and delivery pressure of the compressor. 

b. Pressure at the start of engine compression. 

c. An initial pressure in the bounce space. 

3. Amount and heating value of fuel. 

4. Thermodynamio parameters such as polytropic exponents, 
and ratio of specific heats. 

All of the above parameters may be varied to suit almost any con- 
dition compatible with the assumed geometry. For instance, part load 
behavior may be simulated by merely changing the appropriate pressure 
initial conditions, and then determining the range of fuel energies for 
which the solution is stable. Of course, a different amount of fuel 
applied to the same pressure conditions will lead to a different dis- 
placement-frequency solution. 

The usefulness of the analog computer approach to this problem is 

ihatj 

1. The parameters are introduced as various constants and 
initial conditions and are not manipulated at all, 

2. The oomputer solves the differential equation of motion. 

3. Other quantities may be determined by a minimum amount 
of calculation. 
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It is to be re-emphasized that the analog computer eliminates the 
trial-and-error elements from this problem if the geometry is used as a 
starting point. 
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4. Geometrio Parameters 



There are eight primary geometric parameters and variables asso- 
ciated with the analog approach to the free piston problem. These arei 

1. Area of engine piston, Ag. 

2. Area of compressor piston, A c . 

3. Area of bounce piston, A^. 

4. Engine piston displacement from centerline, x» 

5. Compressor piston displacement from head, y. 

6. Bounoe piston displacement from head, z. 

7. Engine piston position at port closure, Xp. 

8. Piston mass, M. 

Referring to Fig. 1(b), it may be readily seen that 



A b 


“ a e + V 


4.1 


y 


• X + 0, 


4.2 


z 


K 

1 

II 


4.3 



where o end b are the respective clearances of the compressor and bounoe 
cylinders when the engine pistons are touohing. Since the piston is in 
one piece, it is obvious that 

dx _ dy . dz ' 4.4 

dt ~ "31! ” "cTE * 

These relationships are necessary in that x, y, and z and their time 
derivatives are used in the generation of the appropriate forces. 

The mass, M, may be taken as any reasonable value compatible with 
piston size. A computer solution indicating a piston density greater 
than that of steel must be rejected. 
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The treatment of the free piston generator will encompass one-half 
of the system, which will be the right hand piston as illustrated in 
Fig. 3. Positive forces will be considered to be acting to the right. 
Since one of the computer processes will be the summation of these forces, 
the associated sign inversion requires that the negatives of the engine 
and compressor forces be generated. Sinoe one-half of the system is 
being treated, all extensive quantities must be multiplied by two. 



l 
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5o Thermodynamic Parameters. 

Inasmuch as the computer is solving a mechanical problem, all thermo- 
dynamic quantities will be assumed to be given in compatable mechanical 
units. The use of Joule's constant will always be implied whenever 
necessary. Energy leaving a system as heat will be considered negative. 
5ol Pressures. 

The pressure taken as the fundamental variable will be the turbine 
receiver pressure, P^. This will be assumed to remain cyclically constant. 
Using this as a starting point, it is possible to work backwards through 
the system to establish all of the other pressures. These are 



and 







5.1.1 


. & , 


P t 


5.1.2 


n 


he n 3 




- A 


Pc 


5.1.3 


to. 


r e to ri 5 




- P 0 n j 




5.1.4 



p » ii 
r P 

Ps 

P<L 



Pi 



where P p is the engine pressure during scavenging, 

P 3 is the pressure in the scavenge cylinder, 

P^ is the compressor delivery pressure, 

P^ is the compressor intake pressure, 

P 0 is the ambient pressure, and 

r is the ratio of downstream to upstream pressure for 

"e" - Engine exhaust ports, 

"a" - Engine scavenge ports, 

"d M - Compressor discharge valves, and 
"i" - Compressor intake valves. 



5.2 Volumes. j 

In all cylinders, the volume is the product of total displacement 
and area. 

5.3 Temperatures. 1 

Temperature will not be used as such, but rather indirectly in the 

* * 

/ 

determination of the mass of air charged to the engine and the power 



output of the machine. 
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5.4 Mass of Engine Charge. 

The mass of air charged to the engine is not at all a neoessary para- 
meter insofar as the computer is concerned, but must be known in order to 
determine the air-fuel ratio. 

Using the assumption that there is no heat gained or lost in the 
scavenging space, it can be postulated that, upon completion of the sca- 
venging process, ihe engine charge temperature, Tp, will be essentially 
that temperature of the compressor discharge. The mass of engine charge, 
ra e , will then be 

m^ s Pp Vp^RTp, 5«4 0 1 

where Vp is the engine volume at the point of port closure. 

Applying the compression ratio of the compressor to the above equation 
yields 



me s 




ov-l J 



5.4.2 



where T Q is the ambient temperature 

n c is the polytropio exponent of compression in the compressor. 



This result may be expressed in terms of P-fc by applying relationships 
5.1.1 and 5.3.6i 

( 2 * 70 / 7ic~l ) 

f P ^ r A 

o-rr, ' 6 * 4 ° 3 

HX \r e 

If all r 1 s are assumed the same, the above equation reduces to 

V p 






= Yt a U £ sot) 



— 



ftr D 



p^-D ^(+x € -s) 



6 .4 0 4 



Assuming that r is 0.95, and further postulating that n Q is 1.50, the 
above relationship becomes 









0.0113 P t °'\ 



6.4.5 
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where P-^ is expressed in pounds per square inch. For convenience, 
equation 6.4*5 is presented graphically in Fig. 8. 

5.5 Powor Output of the System. 

To determine the power delivered to the turbine, an overall energy 
balance is necessary! 

m d h o -f- ®f HV - (mf + m^h-fc -f Qrej# 5.5.1 

where m^ is the mass of air delivered by the compressor per cycle, 
mf is the mass of fuel injected per cycle. 
h 0 is the enthalpy of entering air. 
ht is the enthalpy of exit gases. 

HV is the lower heating value of fuel ' r /-' 

Qrej is the energy as. heat leaving the system via the cooling 
water per oycle. 

Equation 6.5.1 makes -the customary assumption that the sensible enthalpy 
of the fuel may be negleoted. 

The quantity m^ may be readily determined by noting from the com- 
puter solution the length of the intake portion of the compressor stroke, 
and applying the ideal gas law 

m d = PiA 0 Yj/ RT o# ’ 5.5.2 

where Y^ is the length of the intake portion of the compressor stroke. 

The quantity mf H7 is a known operating variable and Iiq is tabulated. 
The quantity Q re j is composed ofi 

1. Energy as heat leaving the engine. 

2. Energy as heat leaving the compressor. 

t 

3. Energy lost through frictional heating. It will be assumed 
that this will be totally lost to the ooolant. 

The determination of Q re j will be deferred until more speoifio 
differential equations are formulated. ' \ 

Assuming for the moment that Q^j is known, the only remaining un- 
known in equation 5.5.1 is h-fc, the enthalpy of the exit gasses. From 

V' 
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Figure 8 0 Engine Charge Mae* per Cubic Foot Pistcn Displacenent 0 
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this quantity, the temperature of the exit gas may be determined, and 
also the ideal isentropic work available per pound of fluid. The pro- 
duct of this work, the total mass delivered per cycle, and the frequency 
will then give the air horsepower of the system. It is to be noted again 
that this treatment involves one-half of the machine, so that the " total” 
mass mentioned above is twice that given in equation 5.6.2, 



'l 
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6. Electronic Analog of an Internal Combustion Engine. 

6.1 General. 

The force exerted upon the engine piston may be most appropriately 
represented by referring to a typical indicator diagram such as depicted 
in Fig. 4(a) on page 9. Since the free piston device has no crank mech- 
anism, the only fixed geometrical point is that of port closure. The 
thermodynamic state at this point will be known if it is assumed that 
successive cycles are identical. All other point of force as a function 
of displacement are determined by the dynamics of the system. In order 
to simulate this situation it is necessary to obtain a computer solution 
of the differential equation of the first law of thermodynamics. 

In the following analysis, the approximate air-standard is used, which 
assumes : 



1. Air is the working fluid, and 

2 0 The mass of fluid within the cylinder remains constant. 

The latter assumption is also derived from the following two 
practical considerations : 

1. The mathematical complication involved in a changing 
mass system, while possible to solve on 'the computer, would require 
more components than are available. 

2. The mass change is within the same order of magnitude 
as the accuracy of the computer, this there would be no observable 
gain from such a refinement. 

6.2 The General Differential Equation for the Internal Combustion Engine. 

The first law of thermodynamics may be written in differential 
form as 



dQ - d W ■+ dE 5 



6 . 2.1 
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whore dQ is the net energy entering or leaving the system as heat. 
dW is the work done by the system. 
dE is the change in internal energy of the system. 

By taking the time derivative of this equation, and by replaoing dW and 

dE by their respective equivalents in terms of pressure and volume, the 

following equation is obtained. 




By solving for the highest derivative of PJ 



dl 

dt 




(K-l) 



dt 




6.2.3 



Now, the volume is given by 

V - A f x , 

d / _ a cW' 
dt " A£ dt 



Ae • 



6.2.4 



6.2.6 



When equations 6.2.4 and 6.2.5 are substituted into equation 
6.2.2, the final result is 



d& J_ ILzJL 

dt % Ag dt 




6.2.6 



where Pg is the engine pressure. 

This equation may be solved on the BEAC in the following steps. 

1. Integrate dPg to form -Pg with the appropriate initial 

TFT 

oondition, Pp. 



2 . 

3. 



Multiply -Pg by kx to form -kPgx. 



Add k-1 dQ 

Ag ' Zt 



to form 










4. Divide this sum by x. This result is, by equation 6.2.6, 



dPg, which is then fed back to Step 1 

ZT 
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The above steps are represented schematically in Fig» 9, This dia- 
gram also shows considerable computer detail which is explained in the 
following notes. 

1. The inputs to Amplifier 5 are to separate the net dQ into 

cTE 

’’fuel" and "non-fuel' 1 terras in order that these quantities may be varied 

independently. In addition, the "fuel" function is non-zero for only a 

portion of the cycle. The output of this amplifier is dQ net. 

"cR; 

2. The relay in the feedback of Amplifier 1 establishes the 
initial condition, Pp, whenever x is greater than Xp. 

3. The combination of the electronic multiplier and Ampli- 
fier 2 is pre-set to give one-fiftieth of the produot. 

4. The multiplier in the feedback of Amplifier 4 has the 
effect of dividing by x, where-x must always be negative for stability. 
The output of this amplifier is pre-set to give ten times the quotient. 

5. Amplifier 3 is set to give ten times the sum of - 



and 




e <2t, 



in order to compensate for the constants inherent in 
the multiplication and division. 

To gain a better insight into the computer solution of equation 
6.2.6 a few simple qualitative examples will be disoussed in the follow- 
ing three sections. 

6.3 Computer Simulation of a Reversible Adiabatic Process. 

The reversible adiabatic prooess is characterized by dQ equal to 

dt 

zero, in which case equation 6.2.6 reduces to 



dP £ 



at 



£ s — 



fr Pei. 



6.3.1 



It may be easily verified that Pgx k = constant is the solution to the 
above equation, which defines the familiar isentropic process. The 
oomputer solution to equation 6.3.1 appears in Fig. 10 with various 
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Figure 9o Electronic Analog of an Internal Combustion Engine 



Pressure, volts (Arbitrary Scale) 



«* * 




Figure 10o Computer Simulation of an Ison tropic Process 0 
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initial conditions. The input "x" was a displaced sinusoid simulating a 
piston motion at about one cycle per seoond. This was produced from the 
arrangement shown in Fig, 11, where "f is an auxiliary variable solving 



which will be used in the forthcoming discussions of the compressor and 
bounce cylinders. This is an artificial stroke input and will be by- 
passed when the entire free piston analog is assembled, 

The solutions depicted in Fig, 10 can be seen to be excellent for 
all but the lowest values of initial condition. The plotted points are 
the numerical values of the indicated equations, and the solid lines are 
the computer solutions. In all cases, the integration was stopped when- 
ever "x" was greater than 30 volts. The departure from the correct 
solution for the low value of initial condition is due to the fact that 
the inherent electronic errors are magnified whenever the computer is 
working with small voltages. This condition must be avoided if at all 
possible. 

It is to be pointed out that the multiplication of -Pg by x 
involves two quadrants. That is to say, -Pg is always negative, while 
x is both positive and negative. Care must be taken to insure that 
this multiplication yields the same result in both quadrants, A slight 
maladjustment in this multiplier will not return the solution to the 
starting point since dPg in one direction will be different from dPc. in 

aT <nr 

the other, 

6,4 Computer Simulation of a Reversible Polytropio Process, 

The reversible polytropio prooess may be achieved by letting 



the equation ^ - O* Also shown are the rt y" and "z” variables 




6*4.1 
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Figure 11, Generation of Artificial Stroke 



6.4.2 



Upon substitution, aquation 6.2.6 beoomss 

dfc _ _ n Pjjc 

cit x 



It may be easily verified that Pgx 11 - constant is the solution to this 

equation* It should be pointed out that while mathematically dQ is as 

defined in equation 6,4,1, the same result may be achieved on the computer 

by letting _dQ be zero, and applying the coefficient n instead of k after 
dt 2 2 ' 

Amplifier 2 , as suggested by equation 6*4. 2* 

Figure 12 illustrates the computer solution of a few samples of this 
case, again with excellent agreement. At this point, it is important to 
realize that there is no net loss of energy from the system over a cycle* 
The energy lost during compression is restored in the same manner during 
the expansion, or 



■ 0 



This may be easily verified on the computer by integrating equation 
6.4.1. Figure 13 shows displacement, pressure, dQ , and (MQUt as 

Cua dt JW 

funotions of time. As expeoted, spiSk# is zero at -the end of a cycle. 




iT 

6.5 Computer Simulation of a Non-polytropic Process. 

Consider a process which is characterized by a net los3 (or gain) 

of energy. This may be accomplished by letting dQ be some funotion other 

ITS 

than that defined by equation 6.4.1. Suppose, for example, that clg 

dt 

were a small negative constant. The resulting differential equation is 

not easy to solve. However, it can be deduoed that a cyclic change in 

volume will not return the pressure to the starting point, sinoe dPE 

1 dt~ 

is different for each direction of motion. A few computer solutions of 
processes of this type are shown in Fig. 14. With the exception of the 
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Displacement, volte (Arbitrary Scale) 

Figure 12. Computer Simulation of the Reversible 
Polytropie Process „ 
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Figure 13o Rerersible Pclytrcpis Process as a Function of Tin;s 0 
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volta (Arbitrary Seal*) 
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Figur* 14o Cctiputar Simulation of the Nen-Polytropi* 
ProceMo 
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isentropic oase, there is net negative work done. This is, of course, the 
work equivalent of the energy lost by the system as heat# Figure 15 
illustrates the same solutions as a function of time* Here it is noted 
that the net energy, Q, appears as a ramp function. The discontinuity in 
the pressure curve is due to the operation of the initial condition relay 
when the displacement exceeded 30 volts# 

This and the preceding examples serve no purpose in the overall aim 
of this investigation, but are presented in order that the reader may 
beoome familiar with the rather peouliar sort of mathematical manipula- 
tions used in setting up the analog computer. In addition, they may serve 
as useful check points in the event that this problem is set up on a 
different computer. 
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Figure 15 o Non-Polytropic Process as a Function of Time, 
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7 0 Electronic Analog of the Standard Diesel Cyole. 

7.1 General 

Equation 6,2,6, and the associated computer arrangement illustrated 

in Fig, 9 represents the most general oase of an engine oyole, wherein 

any form of heat removal or addition may be simulated by the appropriate 

selection of functions for the dQ terms. The investigation of solutions 

"cTE 

of this type, enoorapassing various methods of heat removal and fuel 
injection, is a very broad subject, and should be treated as a separate 
investigation. It will not be practical to do this within the wider scope 
of the free piston problem. In order that the results from this investi- 
gation may be more readily analysed, a standard diesel cycle will be 
employed for the engine. This has the additional practical advantage 
that excessive computer voltages will be avoided,, 

7.2 The Standard Diesel Cycle, 

The standard diesel cycle consists of the following processes; 

1, Polytropic compression, at exponent n, of the engine charge, 

2, Constant pressure addition of energy over a portion of the 
outstroke, 

3, Polytropio expansion of the engine oharge and combustion 
products, at exponent n' , 

4, Constant volume heat rejection to exhaust, 

5, Constant pressure scavenging process. 

Referring to equations 6,2,6 and 6,4.1, these processes may be described 
mathematically as follows; 

1, Compression process. 
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7.2.2 



2. Combustion prooees. 




and 

(dQ\ _ A g- jL. pi 

Wt/p /r-1 * 7.2.4 



and 



3. 



V*t, 



Expansion process# 

) = - 
y x 









7.2.5 



7.2.6 



The above equations may be readily solved on the computer by applying 
the following modifications to Fig. 9, on page 26. 

1. Elimination of Amplifier 5. 

2. Replacement of by §-n when the velocity is negativej 
and by |n’ when the velocity is positive. 

3. Disconnection of dPg from Amplifier 1 during the pro- 

dt 

soribed interval of oombusion, thus making dPjj zero, 

dt" 

Items 2 and 3 above may be accomplished through the use of relays. It 
should be noted, however, that since the free piston engine clearance is 
variable, the constant pressure combustion proooBS must be initiated by 
a change in sign of velocity and terminated when a prescribed amount of 
energy has been added rather than at particular displacements. The 
reason for this is that a change in clearance should not affeot the 
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amount of fuel energy added. 

Figure 16 illustrates the oomputer arrangement to simulate the 
standard diesel cycle. The functions of the additional amplifiers and 
relays are as follows t 

1, Relay 1, actuated by a ohange in sign of velocity, applies 
the coefficient -|n when velocity is negative, and the coefficient i^when 
the velocity is positive, 

2. Relay 2, also actuated by a ohange in sign of velocity, dis- 
connects dPg from Amplifier 1 when the velooity goes positive! and connects 

diT 

the circuit so that dPg may be reconnected by Relay 3 at the proper time, 

dtT 



3, Amplifier 5 performs the integration of 1_ n* PgX, whenever 

. . . 10 
x is positive. The quantity _yi 'P v may be rewritten as 

10 ® 



7o 



P,i 



10 



K JlK-1 



*]= 



Kzjjda 



7,2.7 



This integration, therefore, yields a result proportional to the total 

instantaneous fuel energy. This is oompared with the desired value of 

Qf in Amplifier 6 whioh actuates Relay 3 when (/dQN <jt equals m^H7. 

jUt/n 

Figure 17 illustrates a typical solution using this arrangement. It 
should be pointed out that the operation of this analog does not depend 
upon a pre-set stroke or clearance. 
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Figure 16 0 Electronic Analog of a Standard Diesel Cycle 



V, 



-t- 







Displacement, volt* (Arbitrary Seal*) 



Figure 17 0 Computer Simulation of a Diesel Cycle 0 
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8« Energy Relationships in the Engine 



In the operation of the analog of the standard diesel engine p it will 
be necessary to know the net energy lost as heat from the system; and the 
net energy added during the combustion process# These quantities may be 
determined by the integration of equations 7#2 # 2 f 7#2#4, and 7#2#6 over 
the appropriate displacements# These aret 

1# Energy as heat transferred during compression 



Q*- wi&i - 1 

2. Net energy as heat transferred during combustion 



V , 



Qp = -ArUVi-v*). 



K - 1 

3. Energy as heat transferred during expansion 

py ‘ 



Qx' ■ ~iZC 



K-h! 






«[ftr - 4. 



8.1 



8.2 



8.3 



r 



where the subscripts 1,2,3, and 4 refer to the points illustrated in 
Fig. 17. It can be seen that with the oomputer furnishing a picture 
similar to Fig. 17, all of the required pressures and volumes will be 
known, thus all of these Q* s may be determined. The sum of Qa and 
will be that portion of Q r# j ohargeable to the engine; and Qp is an 
idealized m^HV product whioh does not include heat lost during combustion, 
nor energy lost through incomplete combustion. 
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9. Electronic Analog of a Reciprocating Compressor 



The force exerted upon the compressor cylinder may be most conven- 
iently represented by referring to a typical indicator diagram for a 
compressor, such as appears in Fig. 4(b), page 9. In the case of the 
free piston device, there is no fixed geometrical point. The only known 
quantities are the pressure limits. 

In a manner similar to that used in the engine equations, the 
differential equation of the first law for the compressor may be written 
as 



where A Q is the compressor piston area. 

P c is the compressor pressure. 

y is the compressor displacement from the head. 
dQ c is the rate of heat transfer. 

dt 

This equation, however, is valid only where the mass is constant, not 
during the delivery or intake processes. For these processes, it can 
be said that 

P 0 z P(j# for delivery, and 



P 0 z Pi, for intake. 

The problem of the changing mass system of the compressor may be 
readily solved on the computer with essentially the same arrangement as 
that usod for the engine, but with the use of diode limiters on the 
amplifier generating P 0 . Suoh an arrangement is shown in Fig, 18, 
which solves the poly tropic equation 



d fj 

eft 




9.2 



* 
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It oan be seen from Fig. 18 that -P c is always negative. With this 
being the case, the diodes in the feedback of Amplifier 7 will function 
as follows* 

1. With y negative (instroke), dP will be positive, and-P 0 , 

dt 

the output of Amplifier 1 , will tend to go more negative. This will 
continue until-Pg overrides the bias on Diode 1 at whioh time this tube 
conducts, thus stopping the integration with _ P 0 held oonstant at-P^. 

This situation will continue as long as y is negative, 

2. When y changes sign, dP will beoome negative and the out- 

dt 

put of this amplifier will beoome less negative. When this happens, 
the initial condition will be -P^, whioh has been imposed by Diode 1. 

3. The absolute value of P c will continue to deorease until 
the bias on Diode 2 overrides -P c , This tube will then conduct, thus 
holding the output at-P^. This will continue until y changes sign and 
the integration is started at the initial condition -P^. From this point 
on the oycle will repeat. 

Thus it oan be seen that by the use of two diode limiters, equation 
7.1 may be solved in the case of a compressor without taking into con- 
sideration the changing mass problem. Furthermore, this arrangement 
will simulate a reciprocating compressor that has a variable stroke, 
whioh is precisely the situation in the free piston maohine. This is 
beoause the integrations covering the compression and expansion strokes 
are initiated by a change in sign of the velooity, not at a particular 
position. 

The analog illustrated in Fig. 18 was constructed using the assump- 
tion that the compression and expansion prooeed at the same polytropic 
exponent. 
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Figure 19 illustrates two solutions of the compressor problem, 
corresponding to two different stroke lengths. Here it is again noted 
that a pre-set stroke and clearance is not necessary. 
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Pressure, rolta (arbitrary scale) 




i 

Figure 19o Computer Sismletio* of a Reciprocating 
Cappre»»cr 0 
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10. Energy Relationships in the Compressor 



In order to solve equation 5.5.1 it is neoessary to know the amount 
of energy removed as heat from the compressor per cyole. Recall that 
applying the coefficient n c instead of k is equivalent to letting 



which produces a reversible polytropio process wherein there is no net 
loss or gain of energy over a oyole. This will be true only in the oase 
of a oonstant mass system. The compressor, however, oompresses one mass 
and expands a different mass. By integrating equation 8.1 over the com- 
pression and expansion strokes, the net energy transferred as heat per 



where m<j is the mass of air delivered by the compressor per cycle. This 
equation neglects the energy transferred during the delivery and intake 
processes. It can be seen that if m d is zero, or if k * n c , there will 
be no heat loss. 

*See Appendix III for derivation. 
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10.1 



oyole is 



'j’rjK-x c) 
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11. Electronic Analog of a Gas Spring 

The force exerted upon the bounoe cylinder is that of a gas spring, 
i.e.jPz^ ~ constant. In this case, equation 6.3.1 may be rewritten as 




where Pb is the pressure on the bounoe piston 
z is the bounce piston displacement. 

On the basis of the preceding analyses, the computer arrangement 
may be immediately depicted as that shown in Fig. 20. The diode in the 
feedback of Amplifier 10 performs tlx? function of an atmospherio relief 
by not permitting Pb to fall below a voltage corresponding to one 
atmosphere, thus simulating the introduction of additional air into 
the bounoe space. A typioal solution from this arrangement is illus- 
trated in Fig. 21, which also illustrates the functioning of the at- 
mospheric relief. 

In the case of the bounce cylinder analog, extreme caution must be 
used in the adjustment of the electronic multiplier forming the Pz 
product, for the reason outlined in section 6.3. This situation is 
aggravated in this case since there may not be an initial condition 
imposed every cycle. That is to say, the solution will be "floating." 

A very slight maladjustment in this multiplier will cause the solution 
to drift rapidly away from the correct value. This malfunctioning is 
illustrated in Fig. 22. This situation must be avoided. 
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Figure 20 o Electronic Analog of a Gas Spring. 
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Figure 21 o Computer Stamlntion of Gas Spring 
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Figure 22 0 Result cf Malfunctioning of the 
Gas Spring Analogo 
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12. Coulomb and VIbcous Friction. 

12.1 The Electronic Analog. 

The force of friction is aBaumed to be composed of the two familiar 
types* Statio friotion, a constant force always opposing the motion, and 
visoous friction, a force proportional to the velocity, again opposing 
the motion. This total foroe is illustrated graphically in Fig. 4(d). 
This may also be expressed mathematically as 



where C is the Coulomb friction force 

v is the viscous friction coefficient. 

The analog to produce a voltage corresponding to this force is 
illustrated below in Fig. 23. 



Figure 23. Electronic Analog of Coulomb and Viscous Friction 
The amplifier has no feedback except for the two diodes and their 
series resistances. TOienever the input is different from zero, even by 
the smallest amount, the simplifier output will tend immediately towards 
saturation. This tendency is checked at the voltage "C" by one or the 
other of the diodes, which when conducting, places the feedback resist- 
ance in the circuit. Figure 24 illustrates the computer analog of the 
force of friction both as a function of velocity and of displacement. 




12.1 



J W.— (^ 1 ) — * . 
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Frictional Force, .volts (Arbitrary Scale) 




Velocity, volts (Arbitrary Scale) 




Displacement, volte (Arbitrary Scale) 



Figure 24o Ccnputer Simulation gf Coulomb ud 
Viscous Friction,, 
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when x is negative. 



12.2 Energy Dissipated by Friotion 0 

The last remaining item as yet unknown in equation 5.5.1 is that 
portion of Qrej chargeable to friotion. 

The force of Friction is given by 

£ - C - vx , 

end 

F f . - -C-vi, 

In order that this force may be integrated dirootly, it shall be assumed 
that the free piston has the approximate motion 

s 



12.2.1 



when x is positive, 



X — C e + -jr/Go^ JrrSt 3 



12.2.2 



where o e is the engine clearance 
f is the frequency 
s is the stroke 
The work done by friotion is 

V14 



or 



'? = F & cbjL , 






Wp= 5 



12.2.3 

12.2.4 



Upon substituting and integrating over a oyole, the net energy lost through 
friotion work is given by 

,3. 



Wr - 9Cs + 

r 



12.2.5 
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13. Computer Operated Relays. 

As outlined in Fig. 16. the use of five relays is required in the 
simulation of a diesel engine. For convenience, these are: 

1 0 A relay to start and stop integration whenever x is greater 
than Xp, the port position. This haB the effect of putting the "tail" on 
the engine picture, thus simulating the scavenging prooess. 

2. Three relays necessary to make dpg zero for a portion of 

dt~ 

the outstroke, thus simulating the oonstant pressure addition of fuel. 

3. One relay to ohange exponent from compression to expansion. 
The above relays, of course, must be actuated by quantities available in 
the computer. In all oases the relays were driven by the output of an 
operational amplifier by use of the scheme illustrated below in Fig. 25. 




Figure 25. Use of Operational Amplifier to Actuate Relays. 

As can be seen above, the amplifier has no feedbaok except diodes to 
limit the output to a voltage, ^sufficient to actuate the relay. The 
action of the above arrangement is to yield an immediate step voltage to 
the relay whenever the algebraic sum of the inputs passes through zero 0 
Figure 26 illustrates the output of suoh an amplifier aotuated by a 
sine wave. It can be seen that this 'output is very nearly a square 
wave. 
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Figure 26* Operation of Computer Operated Relay* 



14. Electronic Analog of the Free Piston Generator. 

The preceding sections may be summarized by referring to an overall 
computer arrangement for the solution of the piston motion. This is 
illustrated in Fig. 27. This diagram is basically the same as Fig. 7, 
except that the details of Figs, 16. 18, 20, 23, and 25 are included. 

The numbers on the amplifiers correspond to those used in the above 
mentioned figures. 



56 




57 



Figure 27 0 Electronic Analog of the Free Piston Gas Generator, 



16* Time and Magnitude Soaling# 

16.1 General, 

All of the preceding seotions have been rather qualitative in nature 
in order that the principles of the analog approach to the problem may be 
presented more clearly. In order to obtain a numerical result, however, 
it is necessary to consider magnitude and time soaling within the oomputer. 

15.2 Soaling the Pressure Analogs. 

As outlined in seotion 6.3, the various pressure analogs were con- 
structed by using as a driving foroe a separately generated displaced sine 
wave, operating at a fixed frequency, based on real time. Such will not 
be the case when each of these analogs becomes a part of the overall 
hook-up as illustrated in Pig. 27. In this case, the integrations form- 
ing displacement and velocity will seek some natural frequency, based on 
Oomputer time, not real time. To this end, a time soaling factor of 100 
will be used. Thus, a real machine operating at 600 oycles per minute 
will appear on the computer as operating at 6 oycles per minute. The 
only criterion necessary in making the pressure integrations compatible 
time-wise with the displacement integrations is that a oommon time soale 
be used. When this is done, all integrations will proceed at the same 
rate. 

The- magnitude soaling is entirely arbitrary, depending upon the 
maximum value of the variable, and its relationship to the 50 volt limi- 
tation of the BEAC amplifiers. The only criterion here is that the same 
scale factor must be applied to the initial conditions* 

There is a definite relationship between -tile various soale factors 
involved within each computing loop simulating the engine, compressor, 
and bounce pressures. The oompressor, for example, solves the equation 



58 




or, in terms of computer voltages 



15.2.1 



O/ » P 










15.2.2 



where o( ' L is ttie magnitude scale, units of n i" per volt, and the barred 
quantities represent the oomputer voltages corresponding to the real 
quantity. By comparing the above equations, it can be seen that O^n 
may not be arbitrarily assigned, but is given by 



oci = 



C/o 



15.2.3 



Similar expressions exist in the engine and bounce analogs. These are 







3 



15.2.4 



and 




or* 



15.2.5 



As outlined in section 6.2, there is no liberty selecting the 
magnitude soales involved in the multiplications and divisions. This 
scaling has already been performed for all amplifiers except those in- 
tegrating to form the various pressures. 

In the oase of the engine, for example, the integrating amplifier 
may be treated as follows i 
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In this amplifier. 







Jt 

Of* 



+ « 



16.2.6 



where is the time scaling factor, oomputer time . 

real time 

When the above constants are associated with the circuit parameters as 
outlined in Appendix H, the following equation is formed 



= a, e 

C ^P e R-i £f/ 

When equation 15.2.3 is substituted, 

P,, _ Q<j_ 

ft, Cjj CX x 



16.2.7 



15.2.8 



By similar arguments, the circuit parameters involved in the com- 
pressor and bounce pressure integrations are given byi 



and 



, 



cXy o( t 



0 



&/0 





C*2 OC t 



15.2.9 



16.2.10 
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The subscripts 1, 7, and 10 used in the above equations refer to the 
number of the amplifier as given in Fig, 27, 

Since x, y, and z refer to the same motion, it must be true that 

o<f x - cSy - c/g , 15.2,11 



and 




15.2.12 



When this substitution is made, there finally results 



<2, ^ Sy 

ft> fty C±y 



‘ 10 



ft/o C&/0 



o( x <X^ 



15.2.13 



This relationship specifies the definite time-magnitude relationship 
between the pressure integrations and the displacement integration. 

15c 3 Scaling the Displacement Integrations,: 

The displacement integrations may best be analysed by redrawing a 
portion of Fig. 27 i 




The following may be written for Amplifier 14. 

* = - jr[-A £ °<r/ £ - Ac°<{>'% + A b c<p b P b + o( Ff 15.3d 
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Prom which 



II 

*• 


A e °^P £ 

/1<x* 5 




A*. 01 * . 

M Oiy ’ 


f?4-i+3 

ff/*3 


Ab °<p B 
M (Xy 3 


and 




R&iw 




Rw 


M 0( x 


The following may be written for Amplifier 1§ 


£ 

x-l 

II 

1 




or 




<2 iS 


or* 


R,sC ils 


Oty. t: 


And for Amplifier l6j 




od x X - - 


(emLjt , 
J <x± J 


or 




£/ 6 


o/y 




C*' X 


By comparing this equation with equation 15.2.13, 



15.3,2 



15.3.3 



16.3.4 



16.3.5 



15.3.6 



15.3.7 



16.3.8 



15.3.9 



seen that 
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♦ 15.3.10 



- <3/o SL f0 

fix*, Hy C & y ft l0 Cj, 0 ~ft /c ^ /c 

This equation means that the RC time constants in the x, Pg, ^ , and 
integrations must be the same. 

* 15.4 Scaling the Friction Analog. 




t 



Again referring to a portion of Fig. 27? 




The equation for this amplifier is 

F & - ±C ~ vo x j 



15.4.1 



from Tdiich 



& 



/3 ^f/3 
/?I3 



V 

<** 



15.4.2 



No further scaling is required in the other amplifiers of Fig. 27. 
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16. The Pescara Machine on the Analog Computer. 

16,1 Results. 

Reference 1 contains the data and a digital analysis of a free piston 
gas generator manufactured by the Pescara organization in Franoe, Appen- 
dix IV contains the applicable original data and some of the results given 
by London. It is interesting to note that the analog approach to the 
problem is nearly the reverse from the digital approaoh. The digital 
approach uses as a starting point an assumed stroke and clearances, and 
then proceeds to calculate some of the geometry, frequency, power output, 
and other quantities of interest. The analog computer approach uses as 
a starting point a particular geometry, and yields as a solution stroke, 
clearance, and frequency. Unfortunately, the data required by the com- 
puter were not entirely oomplete in the original Pescara paper. 
Conspicuously absent were : 

1. The piston mass, M. 

2. The absolute bounoe clearance, b. 

3. The initial bounce pressure^P^, 

All other items were either given or were determined by simple calcula- 
tion. 

Of course, the computer cannot produoe a solution in the absenoe of 
the above items. London, however, had calculated these items in his 
analysis, and his results being the best available data, were used for 
these three parameters. It should be pointed out, however, that this 
digital analysis is not quite compatible with the computer solution in 
that: 

1. The digital analysis used a modified Otto engine cycle, 
while the computer uses a diesel cycle. 

2. The digital method charges friction individually to the 
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engine and compressor in the fora of meohanioal eff ioiencies, while such 
is charged to the machine as a whole in the computer approaoh. 

3. The digital method uses "lumped" quantities suoh aa net 
engine work, net compressor work, and net friotion work, Suoh is not 
possible on the computer sinoe the computer must consider differential 
items in solving the equation of motion. 

With these exceptions noted, the following compares the computer 
solution* of the Pescara maohine with the combined original data and 
London's analysis 

1, Geometrical Results 



Item 


Analog 


Pescara 


Frequency, cyoles/min. 


600 


613 


Stroke, inches 


17.5 


17.6 


Engine effective stroke, inches 


9.74 


9.70 


Engine clearance, inohes 


1.10 


1.14 


Compressor clearance, inches 


2.33 


2.38 


Bounoe Clearance 


2.54 


2.50 


Mass Results 






Compressor air delivered 
lbs/side/cycle 


0.500 


0.380 


Ratio, engine/comp, air 


0.370 


0.616 


Air-fuel ratio 


22.7 


32.0 


Gas flow rate, lbs, /hr. 


36,600 


29,000 


Fuel flow rate, lbs, /hr. 


689 


460 



Refer to Appendix TV for scaling of Pescara Machine 
^Results and Calculations appear in Appendix V. 
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3, Thermodynamic Results 



Item 


Analog 


Pescara 


Delivery Temperature 


1550°R 


1405 °R 


Air Horsepower 


1200 


1138 



The net energy balance wa^ also calculated 0 This showed a discrepancy 

of 6.&/o between the engine work and the compressor plus frictional work 0 
Figures 28, 29 , and 30 are the computer solutions of the engine 
pressure, compressor pressure, and bounce pressure, respectively, as 
functions of their respective displacements, 

16,2 Discussion of Results, 

As shown in the preceding section, the geometrical results v/ere 
in excellent agreement with the actual machine data, 'The rest of the 
results, however, were close, but not as good as might be expected. 

This is attributable to two causes: 

1, The analog computer was actually solving a different kind 
of probl® than the Pescara machine, even though the geometry was essen- 
tially the same. This is because: 

a. A diesel cycle was used instead of a real engine cycle, 

b. Different thermodynamic constants were used, 

c. The frictional constants were assumed at much too low 

a value, 

2. Computer error. In this regard it is to be noted that 
extreme difficulty was encountered in making the bounce analog stable. 

This phenomonon resulted from the error in forming the Pz product by an 
electronic device. Since z changes sign it is imperative that* 

Pz • -P(-z), so that the solution will reproduce itself regardless of 
direction of motion. The electronic multiplier is not capable of the 
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Pressure, psia 



v* 




' Figure 28 0 Genputer Sisruletion of the Bingine of the 
Pescara M»chine 0 
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Pressure, psia Pressure, psia 




Figure 29o CoripMter Simulation of the Conpreseor 
of the Pescara M& chime 0 




Bounce Piston Displacement, inches 



Figure 30 o Computer Simulation of the Bounce 
of the Pescara Mrchine 0 
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accuracy that is required in the case of the bounce piston* The engine 
and compressor analogs are not sensitive to this effect in that an initial 
condition is imposed each cycle which prohibits the solution from drifting 
away, Thu3, the analog of the bounce piston exhibited a phenomenon 
similar to, but not as pronounced as, that shown in Fig, 23, page 51, 

The maximum drift was of the order of 3 volts. The scaling, however 
was such that a 3 volt departure represented a force error of about 
59,000 pounds, which i3 sufficient to disrupt the entire analog, causing 
the solution to become unstable after about three cycles* Thus the 

t 

solution as outlined in Appendix V approaches, but may not be, the steady 
state* This effect may account for most of the discrepancies. 



t 
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17. Conclusions 



The free piston engine problem may be solved on an eleotronio 
analog computer, if such a computer is capable of generating a stable 
solution of the gas spring. The computer can produce a solution to the 
equation of the first lavr of thermodynamics in the simulation of an in- 
ternal combustion engine and of a reciprocating compressor. An analog 
computer is not capable of solving the gas spring problem to the required 
accuracy in this manner due to small inaccuracies in electronic multipli- 
cation# 

Two alternative ways are suggested to produce an accurate analog 
of the gas spring# These ares 

1. Use of a function generator to yield a function which may 
be used to give Pz^ = constant# Such a device would keep positive control 
of the bounce force at all times# 

2# Use of one-quadrant multiplication# It has been pointed 
out that the primary difficulty was that electronic multiplication could 
not give exactly Pz s -P(-z). This difficulty may be avoided by using 
z having always the same sign. This may be accomplished through the use 
of velocity relays and a sign changer to switch both the input and out- 
put of the multiplication such that the multiplier operates only in one 
quadrant# 

It should also be pointed out that the solution of this problem 
required a good sized computer, in that a total of 26 amplifiers, eight 
diode limiters, eight voltage sources, and six electronic multipliers 
were required, in addition to i ‘he necessary read-out equipment# This 
is considered the minimum computer size to handle such a problem# 

Further refinements, or the addition of another compressor analog for 
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the solution of the free piston compressor would, of ooursa, inorease 
the computer requirements. 



y 
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APPENDIX I 



LIST OF ASSUMPTIONS 



1. The scavenge air is furnished at constant pressure and 
temperature o 

2 0 The delivery to the turbine occurs at constant pressure 
and temperature o 

3o The compressor intake is 70°Fj one atmosphere 0 

4° The operating fluid obeys the ideal gas laws# 

5o The pressure drop across each valve and engine port is 

five per cento 

6 0 The geometry is symmetrical with respect the centerline 0 

7o The compressor compresses and expands in accordance with 

a polytropic exponent of 1,30. 

8 0 The sensible enthalpy of the fuel may be neglected. 

9o The engine operates in accordance with an air standard 
diesel cycle, 

10o The energy transferred as heat during the compressor intake 
and delivery is neglected. 

Ho The frictional force consists of coulomb and viscous 
friction,, 

12 0 The friction work is determined in advance by assuming 
that the stroke is sinesoidalo 



Other numerical assumptions regarding the specific problem 
of the Pescara machine are found in Appendix IV. 
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APPENDIX II 



Analog 



Computer Symbols 



(Z) 



Circuit 



Symbol 



Function 



, v, 



J'*” ** 



1* Potentiometer: 
e, 




e A * & e, 



2. Summer: 



<3< Ri 

e.-^rn 

$3 Ri 

e-^AAAAA^ 





3* Summing Integrator: 

C* 



a, ft 

e, — #-A/V\AN 

ft -fAVWW 

<*, <h 

ft 



rlFi 



r 

-1 



>- 







it c>- 



e« 



4, Relay Driver: 
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APPENDIX II (Continued) 



Circuit Symbol Moaning 



5o Boeing Electronic Multiplier* 




Note* Direct sniltiplier output is approximately — fii 

5x 10 8 



6 0 Donner Electronic Multiplier* 





!_ 

IOO 



e,e. 



7 0 Division using Boeing Multipliers 





8 0 Division Using Donner Multiplier* 




Other Special Circuits are explained in texto 
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APPENDIX III 



Ha at Removal from the Compressor 
The compressor cycle is given by the figure below: 




V 



For the processes 1-2 and 3-4, the rate of heat transfer is given by 



dt 



' - t 






\ p cl y 



dt 



It shall be assumed that the delivery and intake processes involve 
no energy transfer as heat» 

For process 1-2 



or 






Upon integrating between the volume limits 1 and Z $ 

a <- iriiitii <■*' r 



i 



Now 



and 



Ac 



i y, ' & ‘ 

p a< [ 1 " ft)”' 7 • 



p, v, -- 7 k, a % . 



Therefore 



$ 



/f- ?!< 






(K-l)l 71c -1) 
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- t). 



By a similar process, Q 34 is gi-ron by 




K ~ he 

(K-i)hc’i) 






The net energy transfer is 









. (LJlz 

[K-Dhc-l 





.5 



now 

7r\ 3 - tk, - - 7hjL, 

Therefore 




— K~ 7} c. 

(K-t)('7>; -1) 





n 

(r-D(*c-i) 
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APPENDIX IV 



PROGRAMMING THE ANALOG COMPUTER FOR THE PESCARA MACHINE 

X* Data. 

Listed below is the combined data from the original data and the 
results given by London;^) 

A. Geometry; 



A e 


0 o 979 ft 2 


(*) 


A c 


5.73 ft 2 


(*) 


A b 


6.77 ft 2 


(*) 


M 


2080 lb 


(*) 


X P 


10.84" 


(t) 


0 


1.24" 


(t) 


b 


21.14" 


(*) 



B. Operating Conditions; 





P* 


4.85 


(*) 




Pt 


64 . 5p8 ia 


(*) 






293 psia 


(*> 



Assumed Quantities; 




C 


50 lb 


▼ 


0.1 lb se c/ft 


n 


1.30 


w n » 


1.40 


i n 0 


1.30 



♦Original Data. 

fCalculated from original data 
4 From London’s solution. 



78 



II. Assignment of Soale Factors. 

A. Pressure Scaling. 

1. The Engine. 

Pp z 64.5 psia 

The compression ratio is estimated as 10 
Therefore Peak pressure 1500 psia 
Let <^-f B z 30 psia/volt = 4320psfa/volt 
Then P^ s 2.15 volts. 

2. The Compressor. 

a 14.0 psia 

Pjj * 14.0 x 4.85 • 67.9 psia 

Let s lo psia/volt s 1440 psf a/volt 

P~j - 1.40 volts 

Pci z 6.79 volts 

3. The Bounoe. 

Let O = 20 psia/volt - 2880 psfa/volt 
Pg » 14.65 volts 

B. Displacement Scaling. 

Let oi^ * per volt * l/24 ft per volt 
Then^xTp » 21.68 volts 
o’ - 2.46 volts 
b = 42.28 volts 

C. Velocity Soaling. 

It is known that 

x dz 1 9 cos 62. 8t inche» 

*max * 50 ft / S90 
Let o(. * 2 ft/aeo/volt 
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D. Acceleration Scaling, 

^Saax 'K 3000 ft/seo^ 

1st Ot-j x 100 ft/soo 2 /Volt 
S, Tim© Soaling. 

let s 100 seconds computer time/l second real time 
F, Friction Force Soaling, 

let <*f F - 20 lbs/volt 



then C 


- 


2.5 volts 


0. Summary. 








= 


4320 psfa/volt 


°*P C 


r 


1440 psfa/volt 




- 


2880 psfa/volt 




- 


l/24 ft/volt 


<*x 


mm 


2 ft/seo/volt 


<*V 


= 


100 ft/seo^/volt 


CXF' 


= 


20 lb/volt 




mm 


100 seo/sec 


Determination 


of 


Circuit Quantitie; 



A, Pressure and Displacement Integrators (Amplifiers 1, 7, 10, and 
16.) 

From equation 15,3. 10» 

RiQii ftioCfio 



& 



/(. 



oC 



Ru 



K Of x CX t 



= 0 . 480 . 



Let R’ s 2 1 meg 
C* s » 1 y< f 
a’s • 0.480 
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v ■<-' 



B. The velocity Integration (Amplifier 15) 



From equation 15.3.7 

&J5 _ 3 10 . — , 0.500 

ft IS ('J-IS c *x °*"t 

Let 

fits "4 *»eg 

a /5 = a ^ 

C. Summation of Foroes (Amplifier 14) 

From equations 15.3.2 to 15.3.5.' 



AW 


Ae °<r c _ 


Q.C5 + 


Hiil 


1 

■ V 

4 

i 




P c 

Ho>: 


±.378 


<2-^3 ftilH 

ftl+3 


A 6 °<p a 
M0Cy 


3.0SL 


e2./W R$-/r 


- °^> 


0.0031 


H/ft 


M <** 


Let: 






~ ^ 
A? /y / ~ 


7)1 ej 

Trupc^ cL/y/ ~ 


0.c,S4' 


= 




O.C39 


ft )+3 ~ ' 1 


Tne* f } iP^. cl /f3 - 


0 . 300 . 


Rio- - to 


j ~tPjsr >- <)-i /y — 


0 , 03 / 0 . 
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D. Formation of Frictional Force (Amplifier 13) 



From equation 15.4.2 



< 2-/3 Hi 



&/3 



R 



13 



<*F h . 



O-Oi. 




R&3 



O- f e j 

10 eg., T &- , 3 - 1-00 
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APPENDIX V 



'» * 



CALCULATION OF RESULTS FROM THE ANALOG COMPUTER SOLUTION 



1. Frequency. 

The measured frequency was six cycles per computer 
minute, corresponding to a real frequency of 600 cycles 
per minute. 

2. Stroke .and clearances. 

Stroke from figure 28 - 17.5 inches. 

Engine effective stroke - 9.74 inches. 

Engine clearance - 1.10 inches. 

Compressor clearance - 2.35 inches. 

Bounce clearance - 2.54 inches. 

3. Compressor delivery. 

Fron equation 5.5.2: mg » A 0 Yj/ RT Q . 

' 1 - 

Y. = 29.9 volts « 14.95 in. 

P* « 14.0 psia 
A* « 5.73 ft 2 

Tq = 530° R 

114 “ (14 ’(^:ltUl6)fil) 4,9S | = °- 500 lb 



4. Engine charge mass. 

From equation 5.4.5: 

m 0 /V p « 0.0113 P t 0 * 7 
* 0.209 lb/ft 3 

V* ( 10.84H0.979 ) * 0.884 ft 3 

P 12 



✓ 

r 
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4. (Continued) 



in = (0.209) (0.884) = 0.185 lb. 

0 

5. Ratio of compressor to engine air. 



6 . 



m e /m d 0.185/0.500 = 0.370. 
Mass of fuel. 

From equation 8.2 : m-pHV » P 2 



<W 



7.-51 volts = 1530 psia ( 

£7 = 3.6 volts - 1.80 in. 

m-HV m 1.4 (1530) (144) (1.8) =115,500 ft-lb 

^ 37 ? (IT) 



Assuming a LHV of 18,200 BTU/lb; 



115.500 ~ 0.00818 lb. 

(18.266) (778) 



7. Air-fuel ratio. 



a-f ration 0.185/0.00818 =22.7 



8. Fuel rate. 



Fuel rate= (0.00818) ( g) ( 600) ( 60) = 589 lb/hr 



9. Gas rate. 

Total mass per stroke per side = 0.508 lb. 

Gas rate » (0.508) (2) (600) (60) « 36,600 lb/hr 



10, Friction work. 



From equation 12.2,5; W^> - 2Cs ^ T7^Tf 

1 S 

W^— ( 2) ( 50 ) ( 17,5 ) ( -rr g )(0ol) ( 17,5 ) 2 

f ( 12 ) 2 (TL2^) 



153 ft-m 



11, Heat lost by compressor. 



Q. - 7h* (^~ y 0 f ( P*) 

(a '-sk\-i) L y 

" [ U.es) 0 - 23 * - 1 J 



From equation 10.2, 
= 5170 ft -lb. 



- 1 



12. Heat lost from engine. 

From equation 8.1, Q) = Jj^Lh^L 

fte _„ 




- 8,000 ft lb 



13. Temperature of exit gasses. 

From equation 51 5.5.1: 

m^ho-h nifHV = (m f -i-m d )h t + Q re j 



m^li- -= (0.500) (778) (126.7) 
ELfHY - 

QrcJ - 



49,200 ft -lb 
115,500 ft -lb 
13,300 ft-lb 



from which 



*4; = 151.400 382 BTU/lb 

T 7 ? /8 ) ( 6 • 508 ) 



From the tables: 
T t « 1550O R# 
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14, Air Horsepower 



From the gas tables, an isontropic expansion 
from the state given in item 13 yields an 
ideal work output of 139 BTU/lb. 



AHP _ ( 139) (776) (0.508) (2) (600) (60) 
33,0000 



1200 



15. Work balance* 

i 

1 £ 

Area under the engine curve « 340 volts . 

W Q =,(340) (4320) (0.979) (1/24) ==60,000 ft-lbs 
Area \mder the compressor curve -= 160 volts 2 . 
W c = (160) (1440) (5.73) (1/24) ® 55,200 ft -lbs. 

W f =* 153 ft-lbs. 



Discrepancy: 4,650 ft-lbs 



if 




'f 
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APPENDIX VI 



THE BOEING ELECTRONIC ANALOG COMPTER 
The Boeing Electronic Analog Computer ( BEAC) consists of an equip- 
ment rack, in ■which is mounted a number of operational amplifiers, 
coefficient potentiometers, limiters, multipliers, and voltage souroes. 

The r^.ck also contains a common power supply for the various items. For 
this investigation, two such raoks were used, A photograph of this equip- 
ment is shown in Fig. 31. Also illustrated is the x-y plotting devioe, 
and the time reoorder. The hook-up to solve a particular problem may be 
made in two wayst 

I ■ 

1. Direotly on the front of the computer, or 

2 0 On a pre-patch board, which is then mounted on the front of 
the machine. All the necessary connections are led in beneath the pre- 
patch panel. 

The booing operational amplifier is a high-quality DC amplifier 
halving an open loop gain of approximately 50,000. This amplifier will 
operate with stability with any feedback ratio up to unity. It .may be 
used for summation, integration, differentiation or more oomplex 

I 

operations. The amplifier remains linear as long as the output is less 
than 50 volts. A glow tube oircuit shorts the output at about 51 volts, 

thus this value is an absolute operating limit 0 

\ 

The general operation of the oomputer is quite stable. This 

!r 

stability may be enhanoed further through the use of an automatic bal- 
ancing system, which restriots the amplifier drift to less than a few 
millivolts per month. 

Detailed descriptions and operating procedures may be found in the 

( 3 ) 

manufacturer's data.' 
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ELECTRONIC MULTIPLIERS 
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